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Abstract 

Consider n independent Goldstein-Kac telegraph processes Xi(t),..., X n (t), n > 

2, t > 0, on the real line R. Each the process Xk(t), k = 1 ,n, describes a stochastic 
motion at constant finite speed q, > 0 of a particle that, at the initial time instant t = 0, 
starts from some initial point x® = Aj.(0) € R and whose evolution is controlled by a 
homogeneous Poisson process Nk(t) of rate Afc > 0. The governing Poisson processes 
Nk(t), k = 1,... ,n, are supposed to be independent as well. Consider the linear form 
of the processes X\(t),, X n (t), n > 2, defined by 

n 

m = y, Ufc Afc (t), 

k =1 

where a*,, k = 1,... ,n, are arbitrary real non-zero constant coefficients. We obtain a 
hyperbolic system of first-order partial differential equations for the joint probability 
densities of the process L(t) and of the directions of motions at arbitrary time t > 0. 

From this system we derive a partial differential equation of order 2 n for the transition 
density of L(t) in the form of a determinant of a block matrix whose elements are the 
differential operators with constant coefficients. The weak convergence of L(t) to a 
homogeneous Wiener process, under Kac’s scaling conditions, is proved. Initial-value 
problems for the transition densities of the sum and difference S±(t) = X\(t) ± X<i(t) 
of two independent telegraph processes with arbitrary parameters, are also posed. 

Keywords: telegraph process, telegraph equation, transition probability density, tem¬ 
pered distributions, linear forms, Hamming metric, hyperbolic partial differential equations, 
determinant of block matrix, initial-value problem, sum and difference of telegraph processes, 
Brownian motion 
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1 Introduction 

The classical Goldstein-Kac telegraph process, first introduced in the works [10] and [14], 
describes the stochastic motion of a particle that moves at constant finite speed on the real 
line M and alternates two possible directions of motion at random Poisson time instants. 
The main properties of this process and its numerous generalizations, as well as some their 
applications, have been studied in a series of works [IH9], [TT1 - H5] . [T5lfl6] . [T91120] . [22H29] . [31]. 
An introduction to the contemporary theory of the telegraph processes and their applications 
in financial modelling can be found in the recently published book m 
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In all of the above mentioned works the main subject of interest was a single particle whose 
evolution develops at some finite speed on R. On the other hand, studying a system of several 
telegraph particles is of a special importance from the point of view of describing various kinds 
of interactions (for the interpretation of different interactions arising in physics, chemistry, 
biology, financial markets, see |181 p. 1173]). However, to the best of the author’s knowledge, 
there exist only a few works where a system of several telegraph processes is considered. A 
closed-form expression for the probability distribution function of the Euclidean distance 
between two independent telegraph processes with arbitrary parameters was derived in [IS]. 
The explicit probability distribution for the sum of two independent telegraph processes with 
the same parameters, both starting from the origin of R, was obtained in d3 

In this article we continue studying the problems devoted to the evolutions of several 
particles moving randomly with finite speed. The subject of our interests is the general 
linear form 

n 

m = E a k £ R, cik 0, 7 1,..., 7z, t E 0? (1*1) 

k=\ 

of n independent telegraph processes X\(t),... ,X n (t), n > 2, t > 0, that, at the initial 
time instant t — 0, simultaneously start from arbitrary initial points x° k £ R, k = 1,... ,n. 
Clearly, linear form (11.11) is a generalization of the sum of two telegraph processes with the 
same parameters studied in tm 

The process L(t) defined by (11.11) is of a great interest since it has a quite natural 
interpretation as a weighted linear combination of n independent stochastic sources of the 
telegraph type. From this point of view, the knowledge of the probabilistic characteristics 
of process L(t) enables to study various important problems, such as attaining a given 
level by the process or being in a fixed interval. Besides purely mathematical interest, this 
interpretation justifies the strong motivation for studying such type of processes. 

The paper is organized as follows. In Section 2 we recall some basic properties of the 
telegraph process that we will substantially be relying on. In Section 3 we describe the 
structure of the distribution of linear form (11.11) and derive a hyperbolic system of first-order 
partial differential equations for the joint densities of process (II.ip and of the directions of 
motion at arbitrary time t > 0. From this system, a hyperbolic partial differential equation 
of order 2 n with constant coefficients for the transition probability density of L(t) is given in 
Section 4 in the form of a determinant having a nice block structure whose elements are the 
commuting first-order differential operators of the system. In Section 5 we prove that, under 
Kac’s scaling conditions, the process L(t) weakly converges to a homogeneous Brownian 
motion on R. Initial-value problems for the sum and difference S ± (t) = X\ (t) ±X 2 (t) of two 
independent telegraph processes with arbitrary parameters, are posed in Section 6. 

2 Some Basic Properties of the Telegraph Process 

The Goldstein-Kac telegraph process describes the stochastic motion of a particle that, at 
the initial time instant t — 0, starts from the origin x = 0 of the real line R and moves with 
some finite constant speed c. The initial direction of the motion (positive or negative) is 
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taken at random with equal probabilities 1/2. The motion is controlled by a homogeneous 
Poisson process of rate A > 0 as follows. At every Poisson event occurence, the particle 
instantaneously takes on the opposite direction and keeps moving with the same speed c 
until the next Poisson event occurs, then it takes on the opposite direction again, and so on. 
This random motion was first studied in the works by Goldstein [10] and Kac [Pf]. 

Let X(t) denote the particle’s position on R at arbitrary time instant t > 0. Since the 
speed c is finite, then at time t > 0 the distribution Pr{A"(f) G dx} is concentrated in the 
finite interval [— ct , ct] which is the support of this distribution. The density f(x,t), x G 
R, t > 0, of the distribution Pr{X(t) G dx} has the structure 

f (x,t) = fs{x,t) + fac(x,t ), 


where f s (x,t ) and f ac (x,t) are the densities of the singular (with respect to the Lebesgue 
measure on the line) and of the absolutely continuous components of the distribution of 
X(t), respectively. 

The singular component is obviously concentrated at two terminal points ±cf of the 
interval [— ct, ct] and corresponds to the case when no one Poisson event occurs until the 
moment t and, therefore, the particle does not change its initial direction. Therefore, the 
probability of being at time t > 0 at the terminal points ±ct is 

Pr (X(t) = ct} = Pr {X(t) = -ct} = ^ e~ xt . (2.1) 


The absolutely continuous component of the distribution of X(t) is concentrated in the 
open interval (—ct, ct) and corresponds to the case when at least one Poisson event occurs 
before time instant t and, therefore, the particle changes its initial direction. The probability 
of this event is 

Pr {X(t) G (-ct, ct)} = 1 - e~ xt . (2.2) 

The principal result by Goldstein m and Kac m states that the density / = f(x, t), x G 
[—ct, ct], t > 0, of the distribution of X(t) satisfies the hyperbolic partial differential equation 


d2 f ■ 2X° f 


-c€f = 0, 


dt 2 ' dt dx 2 ’ ^ ^ 

(which is referred to as the telegraph or damped wave equation) and can be found by solving 
with the initial conditions 


f(x,t)\t =o = 6(x), 


df(x,t) 


dt 


t =o 


0, 


(2.4) 


where h(x) is the Dirac delta-function. This means that the transition density f(x,t ) of 
the process X(t) is the fundamental solution (i.e. the Green’s function) of the telegraph 
equation (12.31) . 

The explicit form of the density f(x,t) is given by the formula (see, for instance, [i22l 
Section 2.5] or [2B, Section 0.4]): 


A t 


f(x, t ) = —— [5(ct — x) + 5(ct + x)] 


+ 


-,—Xt 


2c 


A 


d 


A/ 0 ( - Vc 2 t 2 - x 2 ] + —I Q ( -Vc 2 t 2 


A 




Q(ct — |x|), 


(2.5) 














where @(x) is the Heaviside step function 
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@(x) = 


1, if x > 0, 

0, if' x < 0. 


( 2 . 6 ) 


Here Iq(z) is the modified Bessel function of order zero (that is, the Bessel function with 
imaginary argument) given by the formula 


J o{z) = 


k =0 




2k 


The first term in (12.5jl 




f s (x, t) = — [5(ct - x) + 5(ct + x)] 


(2.7) 


( 2 . 8 ) 


represents the density (as a generalized function) of the singular part of the distribution of 
X(t) concentrated at two terminal points ±cf of the interval [—ct, ct ], while the second term 

in m 


f ac (x,t) = 


Xe 


—X t 


2c 


A/ n 


- Vc 2 t 2 - X 2 + 7-/ 0 

c at 


d 


A 


V cH 2 — x‘ 


Q(ct — |x|), 


(2.9) 


is the density of the absolutely continuous part of the distribution concentrated in the open 
interval (—ct,ct). 

Denote by X x ° (t) the telegraph process starting from arbitrary initial point x° G M. It is 
clear that the transition density of X x °(t) emerges from (I2.5K by the replacement x ^ x — x° 
and has the form 


-,-A t 


f x \x : t ) = —— [S(ct — (x — x 0 )) + 5(ct + (x - x 0 ))] 




+ 


2c 


A 


d 


A/q ( — \Jc 2 t 2 — (x — x 0 ) 2 ) + —/q ( — a/ c 2 t 2 — (x — x 0 ) 2 


dt 


A 


Q(ct — |x — x°|). 
(2 - 10) 

The support of the distribution of X x (t) is the close interval [x° — ct,x° + ct\. The first term 
in fl2T0D 


D —Xt 


ft (x, t ) = — [S(ct - (x - x°)) + S(ct + (x - I 0 ))] 


( 2 . 11 ) 


is the singular part of the density concentrated at the two terminal points x° ± ct of the 
support, while the second term 


= - (* - *°)A 


d 


A 


+7^0 ( - A Jc 2 t 2 - (x - X 0 ) 2 


( 2 . 12 ) 


0(c£ — |x — x°|), 


is the density of the absolutely continuous part of the distribution of A^ 0 (t) concentrated in 
the open interval (x° — ct, x° + ct). 


























5 


The characteristic function of the telegraph process starting from the origin x = 0 with 
density (I2.5[) is given by the formula (see [H2, Section 2.4]): 


H(a, t) = e 


—A t 


+ 


cosh (tV A 2 — c 2 a 2 j + 
cos (ty/c 2 a 2 — A 2 j + 


A 


A 


c 2 a 2 


sinh (tV A 2 - c 2 a 2 j 1 {| q |<a}. 


\Jc 2 a 2 — A 2 


sm 


ty/ c 2 a 2 — A 2 j 


{!“!>!} 


(2.13) 

where lr.\ is the indicator function, «6l, t > 0. Clearly, if the process starts from some 
arbitrary point x° G R and has the transition density (12.101) , then its characteristic function 
H x °(a,t ) expresses through (j2 . 13[) as follows: 

H x °(a,t) = e iax °H(a,t), a 6 1, t > 0. (2.14) 


3 Structure of Distribution and System of Equations 


Let Xy(f),... ,Xn n (t), n > 2, t > 0, be independent Goldstein-Kac telegraph processes on 
the real line M that, at the initial time instant t = 0, simultaneously start from the initial 
points xi, ...,xJgR, respectively. For the sake of simplicity, we omit thereafter the upper 

x° 

indices by identifying X r . k {t) = X k k (t), k — 1,... ,n, bearing in mind, however, the fact that 
the process Xk(t) starts from the initial point x k . Each process X^t), k — 1,. .., n, has some 
constant finite speed Ck > 0 and is controlled by a homogeneous Poisson process N k {t) of rate 
A k > 0, as described above. All these Poisson processes N k (t), k — 1,..., n, are supposed to 
be independent as well. Consider the linear form of the processes AA(t),... ,X n (t), n > 2, 
defined by the equality 

n 

m = E a k X k (t), a k G M, a k ^ 0, k — l,...,n, t> 0, (3-1) 

k =1 

where a k , k — 1,... ,n, are arbitrary real non-zero constant coefficients. 

To describe the structure of the distribution of L(t), consider the following partition of 
the set of indices: 

I + = {ii,, i k } such that a is > 0 for all i s G J + , 1 < s < k, 

I~ = {ii,..., i m } such that a it < 0 for all ii G I~, 1 < l < m, k + m = n. 

The support of the distribution $(x, t) = Pr {L(t) < x} of process L{t) is the close interval 
depending on the coefficients a k , speeds c k and start points x° k and having the form: 


supp L(t) = 


E 

k =1 


a k x° k 


ai ‘ Ci ° - H 

^i s £l+ ii&I- 


k =1 


, a ^ x k + 1 ( Y ai ‘ Ci ° ~ 

i s e/+ iiGl~ 


(3.2)' 

In particular, if all a k > 0, k — 1,... ,n, then the set I~ is empty and, therefore, support 
(13.2[) takes the form 


supp L(t) = 


Y a k( X k ~ c kt), Y ak ( X k + Ckt ) 


k =1 


k =1 


(3.3) 
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At arbitrary time instant t > 0, the distribution <f>(a :,t) contains the singular and ab¬ 
solutely continuous components. The singular part of the distribution corresponds to the 
case, when no one Poisson event (of any Poisson process N k (t), k = 1 occurs by 

time instant t. It is concentrated in the finite point set M s = {qi, ..., q 2 n} C supp L(t) that 
contains 2 n singularity points (each qj is counted according to its multiplicity): 

n n 

Qj = ^2a k x° k + t^2a k i J k c k , j = l,...,2 n , (3.4) 

k =1 k =1 

where i J k = ± 1, k — 1,..., n, are the elements of the ordered sequence crj = ... j = 

1,... ,2 n , of length n. The sign of each i J k , (that is, either +1 or —1), is determined by the 
initial direction (either positive or negative, respectively) taken by the telegraph process 
Xk(t). Emphasize that some qj may coincide in dependence on the particular values of the 
start points x k , coefficients a*, and speeds c k . 

Note that both the terminal points of support (13.21) are singular and, therefore, they 
belong to M s , that is, 

n / 

y ] a k x k ^ ^ ( y i a i s c i a — y " 

k =1 ijGJ - 


Other singular points are the interior points of support (13.21) . It is easy to see that the 
probability of being at arbitrary singularity point qj (taking into account its multiplicity) at 
time instant t is 



Pr {L{t) = qj}= 2n , j = 1, • • •, 2 n , 

(3.5) 

where 




A = X>- 

(3.6) 


k =1 


From (13. 5j) it obviously follows that, for arbitrary t > 0, 

Pr {L{t) e M s } = e~ xt . (3.7) 

If at least one Poisson event occurs by time instant t, then the process L(t) is located in 
the set M ac = supp L(t) — M S1 which is the support of the absolutely continuous part of the 
distribution and the probability of being in this set at time instant t > 0 is: 

Pr{L(f) e M ac } = l - e~ xt . (3.8) 

Define now the two-state direction processes D 1 (t ),..., D n (t), n > 2, t > 0, where 
D k (t), k = 1,... , n, denotes the direction of the telegraph process X k (t) at instant t > 0. 
This means that D k (t) = +1, if at instant t the process X k (t) is developing in the positive 
direction and D k (t) = — 1 otherwise. 

Introduce the joint probability densities of the process L(t) and of the set of directions 
{Eh(£),..., D n (t)} at arbitrary time instant t > 0 by the relation 

U{x,t ) dx = f{n,...,i n }(x,t) dx = Pr{x < L(t) < x + dx, Di(t) = A,... ,D n (t ) = i n }. (3.9) 
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The set of functions (I3.9[i contains 2 n densities indexed by all the ordered sequences of the 
form <r = {4,..., 4} of length n whose elements i k , k — 1,..., n, are either +1 or —1. 

Our first result is given by the following theorem. 


Theorem 1. The joint probability densities (13.91) satisfy the following hyperbolic system 
of 2 n first-order partial differential equations with constant coefficients: 


dU(x,t ) 
dt 


C<T 


9fcr{x,t) 

dx 


Xfa- (x, t) + ^ Afc Uw (x, t ) . 

k =1 


(3.10) 


where 


(X {4; • • • i ik—li iki 4+1 j ■ ■ ■ On/j 

& {4, • • ■, 4-i) 4) 4+i) • • •) 4}) 


f’rr — ^ ' Qfe4Cfe, (3.11) 

fc=l 

and A is given by (13.61) . 

Proof. Let At > 0 be some time increment. Let N k (t, t + At), fc = 1,..., n, denote the 
number of the events of 4th Poisson process N k (t) that have occurred in the time interval 
(t, t + At). Then, according to the total probability formula, we have: 


Pr{L(t + At) < x, Di(t + At) = 4, • • •, D n (t + At) = 4} 

n ( n 

= JJ(1 - A fc At) Pr l Lit) + A UkikCk < x, Di(t) = 4, ■ ■ •, A»(t) = 4 


fc=i 


fc=i 

t+Al 


I 1/ 11/ A 

y A fc At (1 - AjAt) — / Pr {L(t) + a k c k {-i k {r k - t) + 4(t + At 

fc=i 7=i i 


- T fc )) < x, 




L>i(t) = 4,..., 4>fc_i(t) = 4 _ 1 , Dfc(t) = -4) Ac +1 = 4+1, • • •, Ai(t) = 4} dr k + o(At) 

n 

= JJ(1 - AfcAt) Pr {Lit) < x - Co-At, Z4(t) = 4, • • •, D n (t) = 4} 


k =1 
n n 

n (1 - AjAt) / Pr (L(t) < x - a k i k c k {2(t - r k ) + At), 

fc=i i=i •/ 

D\(t) 4, • • •) -Dfc-i(t) 4’— 1 ; L)fc(t) 4) dDk+i 4+ij • • •) D n (t) 4} dr k T o(At). 


The first term on the right-hand side of this expression is related to the case when no one 

n 

Poisson event has occurred in the time interval (t, t + At), that is, if 22 N k (t,t + At) = 0. 

k =1 

The second (integral) term concerns the case when a single Poisson event has occurred in 

n 

this interval, that is, if 22 N k (t,t + At) = 1. Finally, the term o(At) is related to the 

k =1 

case when more that one Poisson events have occurred in the interval (t, t + At), that is, if 

n 

22 N k (t, t + At) > 2 (one can easily check that all such probabilities have the order o(At)). 

k =1 









Since the probability is a continuous function, then, according to the mean-value theorem 
of classical analysis, for any k there exists a time instant r k G (t, t + At), such that 

Pr {L(t + At) < x , Di(t + At) = 4, • • •, D n (t + At) = 4} 

n 

= JJ(1 - A fc At) Pr {L{t) < x - Co-At, Z4(t) Z4(t) = 4} 


fc=l 


n n 


At ^2 Afc n (! - AjAt) Pr {L(t) < x - a k i k c k (2(t - r k ) + At), 


k =1 j=l 

i+k 


D 1 (t)=i i,..., -Dfc-i(t) = 4-i, A,(4 = -4, 4t+i = 4+i, • • •, A»(4 = 4} + o(At). 

In view of the asymptotic formulas 

n n 

(1 - AjAt) = 1 - AAt + o(At), At JJ (1 - AjAt) = At + o(At), 


j=i 


i=i 


the latter relation can be rewritten as follows: 

Pr{L(t + At) < x , Z4(t + At) = 4, • • •, D n (t + At) = 4} 

= Pr (L(t) < x - Co-At, -Di(t) = 4, • • •, _D n (t) = 4} 

- AAt Pr (L(t) < x - Co-At, Di(t) = 4, • • • , D n (t) = 4} 

n 

+ At^ A fc Pr{L(t) < x - a k i k c k (2(t - 4) + At), 

k =1 

D[(t) 4,..., D k —\{t) 4—i, -^4(4 4, -Dfc+i 4+i, • • •, -^4(4 4} ~t~ o(At). 

In terms of densities (13.91) this equality can be represented in the form: 


c—Co-At 


c— Co-At 


/ <r (e,t + At)de= / 4(44 dt, - AAt / 4(44^ 


x-a k i k c k (2(t-T?)+At) 


It n 

At^A fc / 

fc=i 4 


4«(44 d£ + o(At). 


This can be rewritten as follows: 


[4(£4 + At) -4(44K = 


x— Co-At 


4(44 


4(44 <4 


v— Co-At 


x-a k i k c k (2{t-T k )+At) 


AAt f 4(44^ + At^A fe f 

4 fc=i 4 


4w(44 <^ + o(At). 
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Dividing this equality by At, we can represent it in the form: 


(f,t +At) - f<r{€,t)\d£ = -c e 


CfT At 


cc—Co-At 




c—Co-At 


—oo 

x-a k i k c k (2(t-r^)+At) 


A [ fA(,t)d( + Y.\ k j 

k =1 




/<r(£4) df 

o(At) 


At 


Passing now to the limit, as At —> 0, and taking into account that T k —>• t in this case, we 
obtain 

X X 

suitt) ^ = _ c „ux,t)-\ J t)d( + j2 A* f /»«(?,*)#• 

—oo —OO fc—1 _ QO 

Differentiating this equality in x, we finally arrive at (13.1011 . 

Since the principal part of system (13.10|) is strictly hyperbolic, then system (13.101) itself 
is hyperbolic. The theorem is thus completely proved. □ 



Remark 1. Note that system (13.10p consists of 2 n first-order partial differential equations, 
however the equation for each density f a (x,t) contains (besides this function itself) only n 
other densities f^(k){x,t), k = 1 ,...,n. This means that each density f a (x,t) indexed by 
some ordered sequence cr = {H,... ,i n } is expressed in terms of n densities f^ k )(x,t), k = 
1,..., n, whose indices a (k ' > = {fx,..., i k _ i, —i k , i k +i, • • •, i n }, k = 1,..., n, differ from the 
index cr = {z l5 ..., i k _i, i k , i k +i, ■ • ■, i n } hi a single element only. In other words, the equation 
for arbitrary density f a (x,t) in (13.101) with index cr = ... ,i n } links it only with those 

densities whose indices are located from cr at distance 1 in the Hamming metric. 


4 Governing Equation 

Let S n = {(Ti ,..., cr 2 n} denote the ordered set consisting of 2 n sequences, each being of 
length n and having the form cr k = if :> = ±1, j = 1 ,... ,n, k = 

l,...,2 n , n > 2. The order in H n may be arbitrary, but fixed. For our purposes 
it is convenient to choose and fix the lexicographical order of the sequences in E n = 
|crfc = {^,4*^)..., in^}, if ] = ±1, j — 1,..., n, k — 1,..., 2 n |, that is, the order 

< 7 ! = {- 1 ,- 1 ,..., - 1 ,- 1 }, 

cr 2 = { — 1, —1,. . ., —1, +1}, 
cr 3 = { — 1, —1,. . ., +1, —1}, 
cr 4 = {— 1 ,- 1 ,.. .,+!,+!}, 


(x2 n — {+lj +1) • • •, +1, Til- 

Note that this lexicographical order is isomorphic to the binary one by the identification 
— 1 i —y 0 and +1 i —> 1, however, for the sake of visuality, we keep just the lexicographical 
order. 
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Let p(-, •) : H„ x S„ 4 {0,1,..., n} be the Hamming metric. For arbitrary element 
cr k G S n , k — 1,..., 2 n , define a subset C H n by the formula: 


M k = {cr s G H„ : p(cr a , cr fc ) = 1}, k = 1,. . ., T. 

Identifying the notations f k (x,t) = f ak (x,t), c k = c CTfc , k — 1,..., 2”, system (13.101) can be 
represented in the following ordered form: 


df k (x,t ) _ _ df k (x,t ) 
<9f k dx 


A/fcOM)+ ^3 A m fm(x,t ), 

{m : <r m gM fc } 


^ — 1, • • •, 2 n - (4.1) 


The main subject of our interest is the sum of functions (13.91) 

2 n 

POM) = ^2fk(,x,t), (4.2) 

fc=i 


which is the transition probability density of the process L{t) defined by (13.11) . 
Introduce into consideration the column-vector of dimension 2 n 


f = f {x,t) = (f 1 (x,t)j 2 (x,t),... j 2 "-(x,t)y 
and the diagonal (2 n x 2 n )-matrix differential operator 

On = diag{ A t , k= 1 ,.... 2"}, 


(4.3) 


where A k , k — 1,..., 2 n , are the differential operators 

d d 

A k = - + c k —, fc = l,...,2". 

at ax 

Define the scalar (2 n x 2 n )-matrix A„ = ||^ sm ||, s, m = 1,..., 2 n , with the elements 

{ A, if s = m, 

-Afe, if p(cr a , cr m ) = 1 and 4^ 7^ 4”^ s, m = 1,..., 2 n . (4.4) 

0, otherwise, 

In other words, the matrix A n has the following structure. All the diagonal elements are 
equal to A. At the intersection of the s-th row and the m-th column (corresponding to the 
sequences cr s = {44 44 • • • ; 4^} and cr m = {4 m \ ■ ■ ■, in^}, such that the Hamming 

metric between them is 1), the element — X k is located, where k is the position number of 
the non-coinciding elements in these sequences cr s and er m . Note that, since the Hamming 
metric between these sequences is 1, such position number k is unique. All other elements 
of the matrix are zeros. From this definition it follows that each row or column of matrix 
A n contains (n + 1) non-zero elements and 2 n — (n + 1) zeros. The sum of all the elements 
of every row or column of matrix A n is zero in view of the definition of A given by (13.61) . 

In these notations the system (14.11) can be represented in the matrix form 


[®n + A n ] f = 0, 


(4.5) 
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where 0 = (0, 0,..., 0) is the zero vector of dimension 2 n . 

Theorem 2. The transition probability density p(x,t) of the process L(t ) given by (I4.2[) 
satisfies the following hyperbolic partial differential equation of order 2 n with constant coef¬ 
ficients 

{Det [£>„ + A n ]}p(x,t) = 0, (4.6) 

where Det [D n + A n ] is the determinant of the matrix differential operator [£> n + A n ]. 

Proof. The proof immediately emerges by applying the Determinant Theorem [T9l the 
Theorem], [20j Theorem 1] (see also [T71 Theorem 2]) to system (14.5ft . According to this De¬ 
terminant Theorem, in order to extract the governing equation from a system of first-order 
differential equations with commuting differential operators, one only needs to compute the 
determinant of this system whose elements are the commuting differential operators (more¬ 
over, this theorem is also true in the case when the determinant is replaced by a polylinear 
form defined on an arbitrary commutative ring over the field of complex numbers). Since 
the differential operators Ak, k = l,...,2 n , commute with each other, this Determinant 
Theorem is applicable to our case. 

The hyperbolicity of equation (14.61) follows from the hyperbolicity of system (13.101) (or 
(14.5ft ). The theorem is proved. □ 


Remark 2. Derivation of a general analytical formula for the determinant Det [£> n + A n ] 
is a fairly difficult algebraic problem which lies apart of the purposes of this article. Nev¬ 
ertheless, this problem can considerably be simplified, if we notice that, from the form of 
system (13.101) . it follows that matrix £) n + A n has the block structure 


D n + A n 


1 + A-n-i 


E 


n—1 


E n _i \ 

1 + A-n-1 J 


(4.7) 


where the blocks in (14.71) are composed of the following ( 2 n 1 x 2 n ^-matrices: 


®n-i = dia g{^fc> k = 1, 


~)Tl — 1 


}, Si 2 ! = img{A k , k = 2”" 1 + 1,.... 2”}, (4.8) 


the (2 n_1 x 2" _1 )-matrix A n _i is dehned similarly (14.41) (but taking into account its dimen¬ 
sion), and E n _i = — AiE„_i, where E n _i is the unit matrix of dimension ( 2 n ~ l x2 n_1 ). Since 
the matrix E n _x commute with + A n _i] (and, of course, with [S d 2 1 + A n _i] ), then 

applying the well-known Schur’s formulas for the even-order determinants of block matrices 
to (14.71) . we obtain: 


Det [® n + A n ] = Det (S^ + A n _ 1 )(2)^ 1 + K-i) - A \E n . 


( 2 ) 


(4.9) 


Formula (I4.9j) reduces computation of a determinant of dimension (2 n x 2 n ) to the compu¬ 
tation of a determinant of dimension (2 n ~ 1 x 2 n ~ 1 ). It is also useful to note that, in view 
of definition (14.41) . for arbitrary n > 2, the relation D n + A n = [D n + A n ] 2 holds. In other 
words, the matrix T) n + A n coincides with its transposed matrix. 

Clearly, this approach can be extended recurrently to the determinants of lower di¬ 
mensions with final obtaining an explicit (but complicated) formula for the determinant 
Det [£> n + A n ], however this is not our purpose here. 
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Remark 3. In the particular case n — 2, we obtain the (4 x 4)-determinant 


In the case n 


Det [D 3 + A 3 ] 




Ai + A 

—a 2 

i -Ai 

0 



Det [2? 2 

+ A 2 ] — 

—a 2 

A 2 + A 

: 0 

-Ai 


(4.10 

-Aj 

0 

: A 3 + A 

—a 2 







0 

—Ai 

—a 2 

A 4 + A 



the following (8 x 8)-determinant 

emerges: 




Ai + X 

—a 3 

—a 2 

0 

i -Ai 

0 

0 

0 

—a 3 

A 2 + A 

0 

—a 2 

; 0 

-Ai 

0 

0 

—a 2 

0 

A 3 + A 

—a 3 

: 0 

0 

-Ai 

0 

0 

—a 2 

—a 3 

A 4 + A 

: 0 

0 

0 

-Ax 

-Ai 

0 

0 

0 

: A 5 + A 

—A 3 

—a 2 

0 

0 

-A! 

0 

0 

: —A 3 

Aq t a 

0 

—a 2 

0 

0 

-A! 

0 

—a 2 

0 

Aj + A 

— a 3 

0 

0 

0 

-A: 

: 0 

— a 2 

— A 3 

A g + A 


From (14.1011 and (I4.11j) we clearly see the nice block structure of matrix D n + A„, as it was 
noted in (14. 7ft . It is also seen that the diagonal blocks of determinant (14.lip are structurally 
similar to (14.1011 . Such determinants can, therefore, be evaluated by applying the recurrent 
formula (14. 9ft . Note that if we take some other order of the sequences {er^, k = 1,..., 2 n }, 
(this corresponds to some pairwise changes of rows and columns in matrix + A n ), the 
determinant Det [T) n + A n ] keeps the same value. 


Remark 4- To obtain the fundamental solution of partial differential equation (14.6H we 
should solve it with the initial conditions 


p(x,t)\t=o = 5 



d k p(x , t) 


dt k 


= 0, k = 1, 


2 n — 1, 


t =o 


(4.12) 


where <5(-) is the Dirac delta-function. The first condition in (14.1211 expresses the obvious fact 
that, at the initial time moment t = 0, the density of process L(t) is entirely concentrated 

n 

at the start point ^ a^x®- 

k =1 


To pose the initial-value problem for the transition density p(x,t) of process L{t) we 
need to find the respective initial conditions. To do this, we may use the known formulas 
(I2.13P and (12.14ft for the characteristic function of the telegraph process. Since the telegraph 
processes X ^(f), k = l,...,n, are independent, then, in view of (12.1311 and (I2.14I1 . the 
characteristic function of their linear form L{t) is given by the formula: 


H L (a,t) = exp ( —At + ia ^ a k x 0 k j H k (a k a, t), 


k =i 


04 G 


t > 0, 


(4.13) 


k =1 
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where 

H k (S,t) = 

+ 


cosh tJXl-cie + 


Aft 


cos h c kt 2 “A fc + 


VK - 4e 

Aft: 


sinh U Xl - de 


kl ' Vcte^ k 

In particular, setting t = 0 in (14.131) we get the formula 


sin \t\ dP~X-, 


hlSIS^} 




(4.14) 


( U 

H L (a, 0) = exp( ia 


a k x° k 


k= 1 


and its inverting yields the first initial condition in (14.12[) . To obtain other initial conditions, 
we should differentiate (in t ) characteristic function (14. 13[) the respective number of times, 
then inverting (in a) the result of such differentiation and setting then t = 0. 


Remark 5. From the hyperbolicity of equation (14.6[) and initial conditions (14.12ft (more 
precisely, from the first inital condition of (14.121) ) it follows that the fundamental solution 
f(x,t ) of equation (14. 6 1) is a generalized function and, therefore, the differential operator 
Det [£> n + A n ] in (14.61) is treated, for any fixed t > 0, as the differential operator acting in the 
space of generalized functions 5?'. The elements of 5?' are called the tempered distributions. 
Such interpretation becomes more visual if we note that solving the initial-value problem 
fl4.6D -f l4.12D is equivalent to solving the inhomogeneous equation 


{Det. [D„ + A„]} f(x,t ) = S(t) S 



(4.15) 


where the generalized function on the right-hand side of (14.15[) represents the weighted 
sum of the instant point-like sources concentrated, at the initial time moment t — 0, at 

n 

the point ^ a k% k - In such form of writing the initial-value problem (I4.15D , the operator 

k =1 

Det [£> n + A n ] : 5?’ —> 5?' is the differential operator acting from 5?' to itself. From this 
point of view, solving the differential equation (14.6D with initial conditions (I4.12D means 
finding a generalized function f(x, t ) G Xf such that the differential operator Det [X> n + A n ] 

transforms it into the generalized function 5(t)6 ( x — Y2 a k x k ) £ XT”. Since the initial- 

V k= 1 7 

value problem (I4.6D - (14.12D is well-posed (due to the hyperbolicity of equation (I4.6|) ). such 
generalized function f(x,t ) exists and is unique in 5?’ for any hxed t > 0. Therefore, the 
fundamental solution f(x,t) of the linear form L(t) dehned by (13.ID is the Green’s function 
of the initial-value problem (I4.6D - (I4.12D . 

The same concerns the initial-value problem for the transition probability density p(x, t) 
and the respective initial conditions are determined as described in Remark 4 above. Such 
initial-value problem can also be represented in the form of a inhomogeneous partial differ¬ 
ential equation similar to (j4.15|) . but with another generalized function on its right-hand side 
determined by the initial conditions for the transition density p(x,t). 
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5 Limit Theorem 

In this section we examine the limiting behaviour of the linear form L(t) defined by (13.11) 
when the parameters of the telegraph processes tend to infinity in such a way that the 
following Kac’s scaling conditions fulfil: 

c 2 

X k —t +oo, c fc -> +oo, - - 6k, Qk > 0, k — 1,... ,n. (5.1) 

It is well known that, under condition (15.11) . each the telegraph process X k (t) weakly con¬ 
verges to the homogeneous Wiener process W k (t) with zero drift and diffusion coefficient 
a‘f. = g k starting from the initial point x k G M. Therefore, it is quite natural to expect that 
the process L(t) converges to the linear form 

n 

W(t) = Y j akW k (t) 

k =1 

of the Wiener processes W k (t), k — 1,... ,n. In the following theorem we prove this fact. 

Theorem 3. Under Kac’s scaling conditions (15.Tj) the weak convergence L(t ) =>■ W(t) 
takes place, where W (t) is the homogeneous Wiener process with the expectation and diffusion 
coefficient given, respectively, by the formulas: 

n n 

EW (t) = Y a A a w = Y Qka k- ( 5 - 2 ) 

k =1 k =1 

Proof. Consider the characteristic function //^(cpf) of the linear form L(t) given by 

fl4T3D: 

( n \ n 

—Xt + ia'Y^akXkj \\_H k {a k a,f), «6i, t > 0, 

k =1 ' k =1 


71 _ 

where, remind, X = X k and functions H k (a k o>,t ) are given by the formula (see (I4.14D ): 
k =1 


H k {a k a, f) 


cosh ( t 




c 2 k a 2 k a 2 


+ 


Xi 




cfafa 2 


sinh ( t 


\[>f 


c 2 k a 2 k a 2 


1 {Ka|<^} 


+ 


cos t 


c\a\a 2 


A 2 fe 


+ 


A& 


y/4 a W 


-A 2 


sin t 


clafa 2 


A l 


1 


{l“ fcQ l > ^} 


J 


k — 1,..., n. 


From conditions (15.11) it follows that 




\ m— 1 


Q k , if m = 2, 
0, if m > 3, 
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for arbitrary k = 1,... ,n. Then the following asymptotic (under conditions (15.11) 1 formula 
holds: 


t\/Xi-clala* = \ k Ull-^ala 




— Afcf 


1 - i 4 aW - 14 41 ala 2 )'- 441 ( 4 I - 


2 A l ^ 


2-4 \Xl 


2-4-6 \A 


_ \ 4 . ^ C fc 2 2. _ 1 „4 4. _ 1 C fc 6 6 j _ 

— 2 \ k akOL ^ 8 A 3 ak01 ^ 16 A 5 ^ 

~ A k t-^t. 

2 

Therefore, taking into account that, under conditions (15.ip . {X k /c k ) —>• +oo, (c|/A|) —* 0, 
for any k = 1 ,,n, we arrive at the following asymptotic formulas: 




~t 1) 




-1 0, 


Thus, 


cosh (^,\JX 2 k - 44a^j ~ cosh ^A fc f - QkCl ^ a t ] , 


A i 


V x l - c l a W 


sinh (^t\jX 2 k - c 2 k a 2 k a^j ~ sinh ^X k t 


g k a 2 k a 2 


t . 


H k (a k a,t ) ~ cosh A k t - 


g k ala 2 


t I + sinh ( A k t — 


gkala 2 


t = exp A k t - 


g k a 2 k a 2 


and the following asymptotic formula holds: 


]~J H k (a k a, t) ~ exp 


k =1 


Xt - ~ ( Y Qk(1 k ) 


,fc=l 


(under conditions (15■ 1 1) 1. 


Therefore, we finally obtain (under conditions (15.11) 1 the convergence 


lim H L (a,t) — ex p 

Cfc ? Afc — ^H-OO 
( c fc/^ k)—*Qk 


ia Y a k x° k - - ( Y Qkal I a 2 t 


k= 1 


yk= 1 


k = 1, ... ,71, 


and the function on the right-hand side of this limiting relation is the characteristic function 
of the homogeneous Wiener process with the expectation and diffusion coefficient given, 
respectively, by formulas (15.21) . From this convergence of characteristic functions, the weak 
convergence L(t) =>■ W(t) follows. The theorem is proved. □ 


Remark 6. One can prove a more strong result concerning the pointwise convergence of 
the distribution of the linear form L(t) to the distribution of W(t) ( that is, convergence 
at every point of the support of L(t )), but this requires a much more complicated analysis 
based on the theory of Cauchy problems in the space of generalized functions. 
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6 Sum and Difference of Two Telegraph Processes 


In this section we apply the results obtained above for studying the sum and difference 

S ± (t)=X 1 (t)±X 2 (t) (6.1) 

of two independent telegraph processes X\ (t) and X 2 (t). 

The sum of two independent telegraph processes on the real line R, both with the same 
parameters C\ = c 2 = c, Ai = A 2 = A, that simultaneously start from the origin 0 6 R, was 
thoroughly studied in D3| and the explicit probability distribution of this sum was obtained. 
It was also shown that the shifted time derivative of the transition density satisfies the 
telegraph equation with doubled parameters 2c and 2A. A functional relation connecting 
the distributions of the difference of two independent telegraph processes with arbitrary 
parameters and of the Euclidean distance between them, was given in [ 181 Remark 4.4], 

The results of the previous sections enable us to consider the generalizations of these 
models and to study the behaviour of the sum and difference of two independent telegraph 
processes Xi(t), X 2 (t) that, at the initial time moment t — 0, simultaneously start from two 
arbitrary initial points x®, ^61 and are developing with arbitrary constant velocities c± 
and C 2 , respectively. The motions are controlled by two independent Poisson processes of 
arbitrary rates Ai and A 2 , respectively, as described above. 

The coefficients of linear form (16.Ill are a\ = 1, a 2 = 1 for the sum S + (t ) and a\ = 
1, a 2 = —1 for the difference S~(t), respectively. Therefore, according to (13.31) . the supports 
of the distributions of S ± (t ) are the intervals 

supp S ± (t) = [(x? ± x°) ~ (ci + c 2 )t, (x? ± x°) + (ci + c 2 )f]. (6.2) 

The lexicographically-ordered set of sequences in this case is 


err — (—1,—1), <J 2 — (—1, +1), cr 3 — (+1, — 1), <r 4 — (+1,+1), 


and according to (13.41) . the support of the sum S + (t) has, therefore, the following singularity 
points: 


Qi = i x i + x\) - (ci + c 2 )t, 
Q2 = ( x i + x 2) ~ (ci - c 2 )t, 
Qs = ( x i + 4) + (ci - c 2 )t, 
qt = {x°i + ^2) + (ci + c 2 )t, 


(terminal point of the support), 
(interior point of the support), 
(interior point of the support), 
(terminal point of the support). 


(6.3) 


By setting X® = x® = 0, c\ = c 2 = c, Ai = A 2 = A, we arrive at the model studied in |T7] 
with the support supp S + (t) = [—2 ct, 2 ct\. In this case formulas (16.31) produce the three 
singularity points, namely, ±2ct (the terminal points of the support) and 0 (the interior 
point of multiplicity 2). 

Similarly, the support of the difference S~(t) has the following singularity points: 


Qi = fa? - x 2 ) ~ (ci - c 2 )t, 
q 2 = ( x °\ - x l) - (ci + C 2 )t, 
Qs = ( x i - A) + (ci + c 2 )t, 
ql = (x? - x° 2 ) + (ci - c 2 )t, 


(interior point of the support), 
(terminal point of the support), 
(terminal point of the support), 
(interior point of the support). 


(6.4) 
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Note that if both the processes Xf(t) and X 2 (t) start from the same initial point x° { = 
= i° £ 1 and have the same speed c± = c 2 = c, then the support of the difference 
S~(t) takes the form supp S~(t ) = [—2ct, 2 ct) with the three singularity points 0, ±2 ct (the 
interior singularity point 0 has multiplicity 2). We see that in this case difference S~(t) 
has the same support and the same singularity points like the sum S + (t) of two telegraph 
processes with the same speed C\ = c 2 = c that simultaneously start from the origin OgI. 
In view of (13.5k 

Pr {S ± (t)=qf} = j = 1,2, 3,4, 

where A = Ai + A 2 . 

According to (13.1 ID . for the sum S + (t ) the coefficients c f = c+ fc , k = 1,2, 3,4, are: 

c l~ = ~{ c l + C 2 ), c 2 = — ( C 1 ~~ C 2); c 3~ = Cl — C 2 , cl = C\+ C 2 . 


Then operators At, k — 1,2, 3,4, take the form: 




.1 d d 

-M + ici ~ C 2 ) di’ 


d d 

A “-m + ici + C2) &' 


(6.5) 


Similarly, for the difference S (t) the coefficients c k = c CT , k = 1,2, 3,4, are: 
C 1 = — C 2 ), c 2 = — (ci + cl), c 3 = Cl + C 2 , c 4 = C 1 — c 2 , 


and, therefore, the operators A k , k = 1 , 2, 3, 4, become 


A _d d 


._d d 

A 2 ~m~ {ci+ C 2 ) ai' 


d d d d 

A ^m + (C1 + C2) a? ^ = m + (C1 “ C2) &' 


( 6 . 6 ) 


The initial-value problems for the transition densities of processes (16 .1 1) are given by the 
following theorem. 


Theorem 4. The transition probability densities p ± (x,t) of processes (16. IK are the solu¬ 
tions of the initial-value problems 


5 + (Al+A2) 


S 2 + 2 (A, + A 2 )4 - 2 (c 2 + 4 )X _ (A, - A2) 2 

( c ? - + ( A ? - A i) )p ± (*,*) = o. 


P ± (x,t) |t=o = s(x- (x?±4)), 
d 2 p ± (x, t) 


dp ± (x, t ) 


dt 


= 0, 


t =0 


dt 2 

9 3 p ± (a:, £) 


<9t 3 


t=o 


t=o 


= ( c ? + c\) r U - (*! ± xjj)), 

— -2(AiC 2 + A 2 c 2 ) S"(x - (x° A pi J. 


(6.7) 


( 6 . 8 ) 
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where 5"{x ) is the second generalized derivative of Dirac delta-function. 

Since equation (16.71) is hyperbolic, then, for arbitrary t > 0, the solutions p ± (x,t) of 
initial-value problems (I6.7[) - (I6.8[) exist and are unique in the class of generalized functions 

y. 


Proof. To begin with, we establish initial conditions (16.8j) . According to (I4.13p . the 
characteristic functions of the processes S ± (t ) are 

H ± (a, t) = exp(—(Ai + A 2 )£ + ia(x 4 ± x®) H 4 (a, t)H 2 (a, t), a6l, t > 0, 

where functions Hi(a,t), H 2 (a,t) are given by (J4T4]). Differentiating H ± = H ± {a,t) in t, 
after some calculations we obtain: 


H ± (a, £)|t= 0 = exp(ia(a;? i^)), 
d 2 H± 

t =o 


dH± 


dt 2 


= _( c 2 + c 2 ) a 2 e -(,; ± , o )) 


dt 
d 3 H ± 


= 0, 


t =0 


dt 3 


= 2(A 1 c? + A 2 c> 2 e ia(l S ±I "). 


t =o 


Inverting these functions in a yields initial conditions (16.81) . 

Let us now derive the governing equation for the transition density of the sum S + (t). 
To simplify the notations, we identify operators Ak = k = 1,2, 3,4, by omitting the 
upper index, bearing in mind, however, that we deal with the operators A^ presented by 
(16.5[) . Thus, according to Theorem 2, we should evaluate determinant (14.101) with operators 
Ak given by (16.5p . To do this, we apply formula (14.9p to determinant (I4.10p . We have: 


Det [(D 2 + A 2 ] — 


Ai + A 

—a 2 

-Ai 

0 

—A2 

A2 + A 

0 

-Ai 

-Ai 

0 

A 3 + A 

—A2 

0 

-Ai 

—a 2 

A 4 + A 


= Det 



0 N 


A3 + A — A2 

' - A? 

A 2 A 4 T \) VO 1 


( 6 . 9 ) 


(Ai + A)(A 3 + A) — (Af — A2) — A2 (Ai + A 4 + 2 A) 

—A2(A.2 + A 3 + 2 A) (A2 + A) (A4 + A) — (A| — Ao) 


where, remind, A = Ai + A 2 . In view of (16.51) . 

d 


Ai + A 4 = 2 


dt 1 


A 2 + A 3 — 2 


d_ 

dt 


and, therefore, we have: 
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Det pD 2 + A 2 ] 

(Ai + A)(A 3 + A) — (A 2 — A|) —2A 2 (^ + A) 


— 2A 2 ()| + A) (A 2 + A)(A 4 + A) — (A 2 — A|) 

= [(Ai + A)(A 3 + A) — (Aj — A|)] [(A 2 + A)(A 4 + A) — (A 2 — A|)] — 4A 2 + A^j 

— {A\ + A)(Al 2 + A)(A 3 + A) (At + A) — (A 2 — A 2 ) [(A 4 + A)(AL 3 + A) + ( A 2 + A)(Al 4 + A)] 
+ (A? - A 2 ) 2 - 4A 2 (| + . 

( 6 . 10 ) 


According to (16.51) . we have 


(A\ + A)(A 3 + A) 
(A 2 + A)(A 4 + A) 









(4 



+ 2c 2 


d_ 

dx 



Therefore, 

(Ai + A)(A 2 + A)(A 3 + A)(A 4 + A) = 


! + A ) ~ {cl ~ cl) h 


d \ 2 d 2 

_ + A -(e? + c 2 )g? 


-4c? 21 (g + X 2 


dx 2 V dt 


2 „2 


4cic; 


d 4 
dx 4 ’ 


and 


(A 4 + A)(A . 3 + A) + (A. 2 + A) (A 4 + A) — 2 


d_ 

dt 


777 + A ) - (c? - c 2 


)21 
1 " 2 W 


( 6 . 11 ) 


( 6 . 12 ) 


Substituting (16.111) and (16.121) into (16.10ft we obtain 

1 2 


Det [J) 2 + A 2 ] — 


d \ 

m + x 


2(A? - A|) 


a 


(C? + <£) 


9a; 2 


4c 2 c 2 


g> 4 

<9x 4 


l + A 


4 ) 


— + A ) - (c 4 + c 2 


cf_ 

dx 2 


-i 2 


- 4c?c 2 


cAr 2 
<9 4 


+ (Af-A 2 ) 2 -4A 2 (^ + A 


d 


dx 4 


- 2(A 2 + A 2 ) ( ^ + a \ + 2(A 2 - A^)(c 2 - C 2 )- q ^ + (^1 _ ^) 2 
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d 4 
dx 4 




proving equation (16.71) for the transition density p + (x,t) of the sum S + (t). 

Comparing now operators (16. 5 p and (16.6f) . we see that A± = AA% = Af, A$ = 
A±, A^ = A^. Therefore, as is easy to see, determinant (14.101) written for operators 
A^, k = 1,2, 3,4, takes the same value like determinant (16.91) for operators At, k = 1,2, 3,4. 
Thus, equation (16.71) is also valid for the transition density p~(x,t) of the difference S~(t). 
This completes the proof of the theorem. □ 


Remark 7. By setting C\ = C 2 = c and Ai = A 2 = 
equation 



d 2 a 

ae +4X Ft 


Ac 2 


A in (16.71) we arrive at the fourth-order 
d 2 \ 

w) p {x ’ t) = 0 


and this result is weaker in comparison with the third-order equation obtained in HO formula 
(4.3) therein] for the transition density of the sum S + (t) of two independent telegraph 
processes, both with the same parameters (c, A). It is interesting to note that in the product 
of these operators, the second one represents the Goldstein-Kac telegraph operator with 
doubled parameters (2c, 2A). 


Remark 8. The form of equation (16. 7 p enables us to make some interesting probabilistic 
observations. We see that the first term of equation (16.71) in square brackets represents 
a telegraph-type operator (containing also the free term — (Ai — A 2 ) 2 ) which is invariant 
with respect to parameters Ai,A 2 and ci,C 2 - In other words, if we change Ai for A 2 and 
inversely, and/or Ci for C 2 and inversely, the first telegraph-type term of equation (16.71) 
preserves its form, while the second one can vary its interior signs. From equation (16.71) we 
can conclude that the sum and difference 5 ,± (t) of two independent telegraph processes are 
not the telegraph processes, however they still contain some telegraph-type component. To 
show that, we represent equation (16.71) in the following form: 



































21 


d_ 

dt 


+ (Ai + A 2 ) 


d 2 d 

—+ 2(A 1 + A 2 )- 


2 (c? + c 2 ) 


dx 2 


(Ai 


»4 


a 2 

'cfr 2 


<9 <9 2 

( A i -A 2 ) - + ( c 2 -c^ ) ^ + 2 ( A 2 -A 2) 


p ± (x, f) = 0. 

(6.13) 

We see that the first term of (16.131) contains exactly the telegraph operator quite similar to the 
classical Goldstein-Kac operator (12.3[) with the replacements A )->• Ai + A 2 and c 2 )->• 2(cf + c 2 ). 
The second term of (16.13(1 represents the product of two heat operators and this fact implies 
the presence of a Brownian-type component in the processes S ± (t). Notice also that in 
this product the first operator in square brackets is exactly the standard heat operator (for 
Ai > A 2 and c± > c 2 ), while the second one represents an inverse-time heat operator (that 
is, with the inverse time replacement t i—>• —t) containing also the free term 2(A 2 — A 2 ). 


Remark 9. Solving the initial-value problem (16.71) - (16.81) and obtaining the transition 
densities p ± (x,t) of the sum and difference S ± (t) of two independent telegraph processes 
with arbitrary parameters is a fairly difficult problem that will be realized in the framework 
of another project. 
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